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Abstract 



We describe an experimental test of whether particle decay causes wave 
function collapse. The test uses interference between two well separated, but 
coherent, sources of vector mesons. The short-lived mesons decay before their 
wave functions can overlap, so any interference must involve identical final 
states. Unlike previous tests of nonlocality, the interference involves contin- 
uous variables, momentum and position. Interference can only occur if the 
wave function retains amplitudes for all possible decays. The interference can 
be studied through the transverse momentum spectrum of the reconstructed 
mesons. 
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In 1935, Einstein, Podolsky and Rosen (EPR) showed that quantum mechanics required 
that wave functions can be non-local [1]. When a system is observed, the wave function 
collapses from one which contains amplitudes for a host of possible outcomes to smaller 
set of possibilities, in accord with the measurement. This collapse is instantaneous; much 
has been written about its superluminous nature. Most studies of the EPR paradox have 
tested Bell's inequality [2] using spin correlations, usually with photons produced in pairs 
[3]. Experimenters measure the spin correlations using two polarizers with a varying angle 
between them. Bell found that models with non-local wave functions and models with hidden 
variables produced different angular correlation spectra. Previous tests of non-locality used 
discrete variables like 'pseudo-spin' for CP violation, as with studies using the reaction 
$ -> K+K- [4]. 

We describe a very different system that, in contrast to the K L K S system, is sensitive 
to the collapse of continuous variables in a wave function [5]. Short-lived vector mesons 
(VMs) are produced with a fixed phase relationship at two separated sources. Even though 
the mesons do not come from a single source, and, in fact, share no common history [6], 
the system acts as an interferometer. The meson lifetimes are short compared to the source 
separation, so the mesons decay before their wave functions can spatially overlap. 

Any interference between the two sources must involve the decay products. Interference 
is only possible between identical final states. With the large phase space for final states, 
interference can only occur if the wave functions retain amplitudes for all possible decay 
channels and angular distributions long after the decay takes place [7]. We have previously 
calculated the interference pattern [8]. This letter will focus on the effects of the wave 
function collapse and Bells inequality-like tests, and sketch an alternate derivation of the 
interference, to emphasize the symmetries of the system. 

Figure 1 shows electromagnetic VM production in relativistic heavy ion collisions at large 

— * 

impact parameters, b. A photon from the electromagnetic field of one nucleus fluctuates to 
a virtual quark-anti-quark pair which elastically scatters from the other nucleus, emerging 
as a real vector meson [9]. Either nucleus can emit the vector meson. The momentum 
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transfers from the nuclei are similar, and they remain in the ground state, so it is impossible 
to determine which nucleus emitted the photon and which is the target. 

The electromagnetic interaction (photon) has a long range, while the elastic scattering 
has a short range, around 0.6 fm [10], far smaller than the < 7 fm radius of a heavy nucleus 
or the typical impact parameter. So, VM production takes place essentially 'on top of the 
emitting nucleus, and the two nuclei act as a two-source interferometer. 

Electromagnetic VM production is studied at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory, where gold ions collide at center of mass energies up 
to 200 GeV per nucleon. Starting in 2006, the Large Hadron Collider (LHC) at CERN will 
collide lead ions at a center of mass energy of 5.5 TeV per nucleon. 

The cross sections were previously calculated [11] using the Glauber approach [9] with 
the photon spectrum given by the Weizsacker- Williams virtual photon method [12]. The 
calculated photonuclear cross sections agree with data to within 20%. The cross sections 
are large, about 10% of the hadronic cross section at RHIC, rising to 50% at the LHC. The 
corresponding production rates, more than 100 p°/sec at RHIC, rising to 230,000 p°/sec at 
the LHC, are large enough that it will be easy to collect adequate statistics to study wave 
function collapse. Already, the STAR collaboration [13] has observed more than 10,000 p°, 
which should be enough to observe the interference. 

The impact parameters for these interaction are large compared with the nuclear radii, 
Ra- For p and uj production, the median impact parameter (b) is about 40 fm at RHIC, 
rising to 300 fm at the LHC; for the J/ip, (b) rises from 23 fm at RHIC to about 50 fm 
at the LHC. All are much larger than R A 7 fm for heavy ions. It is possible to select 
events with smaller (b), but still with b > 2R A , by choosing events where VM production is 
accompanied by nuclear breakup [14]. 

The (b) are larger than the distance travelled by most VMs before they decay. The VM 
lifetimes r range from 4 x 10~ 24 s for the p up to 7.5 x 10~ 21 s for the J/tp. The VM are 
produced with typical transverse momentum ~ 2H/Ra ~ 60MeV/c; at mid-rapidity, the 
longitudinal momentum is zero, so VM have a median decay distance d = 2hcr/RAMy. 
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Except for the J/ip, d < (b); for the J/ip at the LHC, d « (6). 

The final state wave function from ion source % at a time t can be expressed schematically 

^(t)i = exp (-t/2r) \V > +(1 - exp (-t/2r)) \DP > (1) 

where r is the vector meson lifetime, \V > is the vector meson wave function, and \DP > 
is the final state. For stable particles, r = oo, the decay products drop out, leaving a 
conventional two-source interferometer. 

The interference can be seen by examining the symmetries of the system. The total 
amplitude At for observing the VM with momentum p at position r, and time, t, depends 
on the production amplitude A(p, x, t') and a propagator P(p, x, t', r, t) which transports the 
meson from x', t' to r, t: 

A T (p, f,t) = J A(p, x\ t')P(p, x, t', f, t)dxdt'. (2) 

The production amplitude A(p,r,t) depends on the electromagnetic field, E(x,t'), nuclear 

— # — # 

density p(x,t) and the amplitude f(p, k) for a photon with momentum k to fluctuate to a 
qq pair and scatter from a nucleon, emerging as a vector meson with momentum p: 

A(p, x>, t') = f(p, k)p(x>, t')E(x>, if) (3) 

The electromagnetic field at a distance b from a nucleus is a Lorentz-contracted pulse with 
a width 6/7 where 7 is the Lorentz boost. When 7 ^> 1, the electric and magnetic fields are 
perpendicular and the overall field may be represented as a stream of almost-real photons, 
with energies up to H'j/b [12]. The photon amplitude is proportional to E(x',t'). The scat- 
tering amplitude is obtained from data; only its symmetries are important here. Absorption 
of the nascent p° is neglected, but could be included with an additional position-dependent 
variable, effectively modifying p(x,t'). 

At large distances, the propagator for a VM with energy u = y 1 My + \p\ 2 may be 
modelled with a plane wave. Neglecting, for now, VM decays, 

P(p, x, t', f, t) = e i&(f-*)-»{t-t>)) (4) 

4 



The nuclear density is symmetric around the center of mass (origin), giving it positive 
parity, while the antisymmetric electric field has negative parity): p(x,t') = p(—x,t') and 
E(x,t') = —E(—x,t'). With this, the range of integration in Eq. (2) can be restricted to a 
single nucleus: 

A T (p, f,t)= f dxdt'p{x, t')E(x, t') e *(p-^(*-*')) %yv-i _ ffa -k) e -^ 3 \ (5) 

Jy>0 

— * 

The only differences between the two nuclei are the phases ±ip-x and between f(p, —k) and 
f(p, k). The latter is because the sign of p z reduces the symmetry of the system. Of course, 
interference is only significant when \f(p, k)\ as \ f(p, —k)\ which occurs near p z = 0. 

— * — * 

The equation simplifies by defining x = b/2 + x'. The bulk of the cross section is 

— * — * 

from when the photon couples coherently to the target nucleus, i.e. when k ■ x' h, so 
the exponential phase is constant over the nucleus. Then, the maximum transverse and 
longitudinal momenta are Ti/Ra and ^Ti/Ra- Emitted photons are subject to similar limits 
[15]. Near p z = 0, the photon momentum, and the momentum exchange due to the scattering 
are very similar so it isn't possible to determine which nucleus emitted the photon, and which 
was the scatterer; in fact, at p = 0, the two momentum transfers are equal and opposite. 

The electromagnetic pulse lasts a time, b/cy, which may be slightly longer than one 
photon period (h/cu). This partial temporal incoherence will reduce the overall production 
amplitude. There is a pairwise cancellation between space-time volume elements dxdt' at 
positions x and —x, so the interference is not affected. 

With this, p • x — p • b/2 and 

A T (p, r,t)= f dx'dt'p(x', t')E(x>, t') e *(p-— (*-*')) k)e^ /2 - f(p, -k)e~^ /2 ] . (6) 

Jy>0 

We now introduce a few approximations. The amplitude for production from the first 
nucleus is Ai(p,r,t) = J y>0 dxdt'p(x,t')E(x,t')f(p, k) and we define c to be the ratio of the 
amplitudes for production from the two nuclei: 

= J y>0 dxdt'p(x,t')E(x,t')f(p,-k) 
I y>0 dxdt'p(x,t')E(x,t')f(p,k) ' 
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The ratio f(p, —k)/f(p, k) does not vary significantly over the nucleus, so the single nucleus 
production amplitude factors out of the integral in Eq. (6). The transverse momenta do not 
affect c, and k z = My/Ap z . Then, 



A T (p,r,t) = A x (p,r,t) 



(8) 



The amplitude factorizes into a magnitude and an interference term. The pt dependence of 
A\ (p, r, t) is dominated by the nuclear form factors, with the bulk of the production having 
Pt < 2%/Ra- Most of the uncertainties discussed earlier do not affect the interference term. 
The time and z variation in E(x, t') should be largely independent of k. In the soft Pomeron 
model, the photon to VM coupling increases slowly with k and has an almost constant phase. 
At RHIC, a photon-meson term is also present, but the phase of c still changes only slowly 
with k [9]. 

The interference is clearest when p z =0. Then c = 1 and the approximations introduced 
in defining c disappear. The amplitude is A T {p,r,t) = 2iA sin(p- b/2). and the cross 
section is 

a ~ \A T (p,f,t)\ 2 = 2A 2 [1 - cos(p • b)}. (9) 

This formula applies for stable particle emission, such as bremsstrahlung photon emission 
in e~e~ [16] and pp collisions [17]. 

For short-lived particles, the situation is more interesting. We can express the cross 
section in Eq. 9 as a sum of a coherent (interfering) and an incoherent (non-interfering) 
term 

2At[l-(l-r))-cos(p-b)}, (10) 



where 77 measures the degree of decoherence, as has been done for the system $ — > K°K 
[4,18]. Complete quantum mechanical coherence between the two sources here means r\ = 0. 

If the interference is restricted to the time the system spends in the vector meson (parent) 
state, then one would expect partial decoherence. According to Eq. 1, only a fraction 
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exp(— (Myb)/ (out)) of the vector mesons will have survived long enough for the amplitude to 
propagate the distance b between the nuclei before the decay. This scenario thus corresponds 
to 

r] = l-exp( ). (11) 

UJT 

For the p, u, and 0, which all have cr <C< b >, rj m and one expects almost complete 
decoherence in this scenario. For the J/^, on the other hand, cr ^< b >, so the decoherence 
would be only partial. This scenario could therefore be distinguished by observing the J '/\& , 
while for the lighter vector mesons it would essentially be indistinguishable from a scenario 
with no interference. 

However, there is a broader issue. A distant observer sees the decay products from the 
original meson, and most VM have many decay modes. The final state may be written 
\DP >= T<jaj\DPj >, where ctj is the amplitude for final state \DPj >. For example, 
the J/ip can decay to e + e~ or 7r + 7r~7r°, among many possibilities. The probability for any 
specific final state is small, less than 7%. 

If the wave function contains amplitudes for every possible decay mode, at times t ^> r, 
the individual particles may be modelled as plane waves and the wave function is 

tt(t) = jT|^ + ce ip-.(£ d -b/2)j e -t d /2r-iu,t d js.y^r^. > ( 12 ) 

where the decays occur at time and displacement Xd = (p/M v )td from the production 
points. Here, Brj are the branching ratios to different final states, and fyj is the k— particle 
final state 

^ = ^W^-^-^it-t,)]^ > 5{ ^ kUk _ u)5 ^ k _ £). (13) 

Here Xk, oju and pk are the particle positions, energy, and momenta, and jk > includes the 
particle types and angular distribution. The 5 functions impose 4-momentum conservation. 

If the decay occurs before the two amplitudes overlap, then either the wave function must 
retain amplitudes for all final states, or the two decays will produce different final states, 



and the interference must be small. Also, as Eq. (9) shows, amplitudes for different decay- 
times must also be included. Thus, the presence or absence of interference tests whether 
particle decay collapses the wave function. 

The two source terms in Eq. (12) entangle the final state wave functions. The phases 
differ by exp [i(p • b + 5)} where S is the phase of c. Any measurement on one decay product 
at least partially collapse the wave function of the others [3]. Detectors could accurately 
measure either the position or momentum of the k final state particles. 'Accurately', is 
compared to the relevant distance (6) or momentum (Ti/b) scales. By these metrics, current 
and planned experiments measure momentum accurately, but not position. 

The interference pattern, Eq. (9) can be seen in the reconstructed VM p T [8]. Because 
b is not generally measurable, a must be integrated over all b. Figure 2 shows the expected 
Pt spectrum for p° production at mid-rapidity at RHIC [8]. The large dip for p T < h/{b) is 
a distinctive experimental signature. 

Alternately, at least in a gedanken experiment, position sensitive detectors could be used 
to localize the decay to a single nucleus, provided the ion trajectories are known. The decay 
J/tp — > e + e~ produces two relatistic electrons that are back-to-back in the transverse plane. 
For p T = 0, a line between the two measured electron positions will intersect one of the ion 
trajectories. When the electron f>k are perpendicular to b then it is possible to determine 
which nucleus emitted the VM. The nonzero meson p T introduces some uncertainty, but not 
enough to encompass both ion trajectories. For detectors 500 fm from the collision point, 
the pointing uncertainty is 16 fm, less than the (b) ~ 50 fm for J/ip at the LHC. 

As with existing tests of Bell's inequality, two detectors, on opposite sides of the collision 
region could randomly measure either position or momentum. All single-detector measure- 
ments are insensitive to what happens in the other detector. However, when both detectors 
measured position, the production point could be determined, while when both measured 
momentum, a null at pr = would be seen, showing the interference. These two possi- 
bilities are only compatible if the wave function collapses when a measurement is made, 
and not earlier, when the meson decays. For J/ip — > e + e~, the collapse would have to be 
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superluminal. 

One possible source of decoherence is the decay timing. The VM are produced nearly at 
rest, but the decay products may be relativistic. The maximum flight time difference from 
the two sources to a detector is b/c. If the detectors could resolve this time difference, this 
could partially localize the production, reducing the coherence. 

Since the probability of producing a VM in grazing collisions {b = 2Ra) is high, about 1% 
at RHIC and 3% at the LHC [11], multiple VM production is also observable. Multi-meson 
final states will exhibit more complicated entanglements, with possibly new behavior. 

We have described a 2-source interferometer for short lived particles, and showed that its 
description requires a non-local wave function. The observation of interference will clearly 
demonstrate that, after a decay, a systems' wave function includes amplitudes for all possible 
decay modes and angular distributions, and does not collapse to a specific decay mode 
until the wave function is externally observed. Measurements of this interference should be 
available soon. The STAR detector has observed exclusive p° production in gold collisions 
at RHIC [13]. Current data should provide higher statistics and an accurate p° px spectrum, 
probing the EPR paradox for continuous variables. 

This work was supported by the U.S. Department of Energy under Contract No. DE- 
AC-03076SF00098 and by the Swedish Research Council (VR). 
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FIG. 1. Diagram showing ultra-peripheral p° production and decay in heavy ion collisions. 
The nuclear momenta follow the z axis, and come closest at z = 0, when their separation (impact 
parameter), b follows the y axis. 
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FIG. 2. Perpendicular momentum spectra for p° production at RHIC, at p z = 0, for gold on 
gold collisions at a center of mass energy of 200 GeV per nucleon. Plotted are dN/dpx, with 
and without interference. The curves are normalized to 1 for pt = and no interference. The 
calculation assumes that the impact parameter is not measured, so the interference is washed out, 
except for pt < 25 MeV/c. 
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